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Abstract: Theoretical predictions show that at low values of Bjorken x the spin structure function, 
gi is influenced by large logarithmic corrections, ln 2 (l/x), which may be predominant in this region. 
These corrections are also partially contained in the NLO part of the standard DGLAP evolution. 
. Here we calculate the non-singlet component of the nucleon structure function, g^ s — g\ — g", and 

its first moment, using a unified evolution equation. This equation incorporates the terms describing 
the NLO DGLAP evolution and the terms contributing to the ln 2 (l/x) rcsummation. In order to 
avoid double counting in the overlapping regions of the phase-space, a unique way of including the 
Oh! NLO terms into the unified evolution equation is proposed. The scheme-independent results obtained 

Oh, from this unified evolution are compared to the NLO fit to experimental data, GRSV'2000. Analysis 

of the first moments of g± shows that the unified evolution including the ln 2 (l/x) resummation goes 
beyond the NLO DGLAP analysis. Corrections generated by double logarithms at low x influence the 

• i-H i 

■ Q 2 -dependence of the first moments strongly. 

.«?: 

1. Introduction 

During the last years the interest in physics of polarized nucleon has increased sig- 
nificantly Experimental data describing the structure of the nucleon are now available 
in various ranges of the momentum transfer, Q 2 , and the Bjorken variable, x. New 
experiments emerge, and new data are expected to be available in the next future. The 
data from the region of low x will be of special importance, since they may improve 
the estimation of the quark contribution to the total nucleon spin. 

The previous data obtained by the EMC collaboration showed that the total 
participation of quarks in the proton spin was very small. This contradicted theo- 
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retical predictions, obtained from the Ellis- Jaffe sum rule. That sum rule expressed 
the moments of quark distributions in terms of the nucleon axial coupling constants. 
Following the Ellis- Jaffe sum rule, the quarks should participate in about one-fifth of 
the total spin of nucleon. The discrepancy between the theoretical expectations and 
the experimental data has been often referred as the "proton spin crisis". 

There have been several explanations proposed to clarify this. One of them refers 
to the not-yet-measured domain of very low values of Bjorken x, x < 1CT 3 . Theoretical 
predictions show that at low x the structure function of the polarized nucleon, g\(x, Q 2 ), 
is influenced by large logarithmic corrections, ln 2 (l/x), 0,0]. As a consequence, large 
contributions to the moments of the structure functions from this region are expected. 

Analysis of the nucleon structure function at low x, including the resummation 
of the logarithmic corrections, was performed in 0,0,3,^,0, 0,11- The logarithmic 
corrections were introduced for the unintegrated parton distributions in Refs. pl,|9|.[10|. 



II]. Those corrections originated from the ladder and the non-ladder bremsstrahlung 



diagrams B,^|, F2|[l3| . 

The recursive equations for resummation of the logarithmic corrections were for- 
mulated. Afterwards, these evolution equations were completed with the DGLAP evo- 
lution terms calculated at the LO accuracy |§||]. Including the DGLAP evolution was 
neccessary for accurate description of the structure functions in the region of moderate 
and large values of x. 

Since the phase-space domains for the logarithmic and the DGLAP contributions 
overlapped partially, subtractions from the evolution kernels were made in order to 
avoid double counting in those regions. It was shown that including the ln 2 (l/x) 
terms into the evolution equations generated more singular behaviour of g\(x,Q 2 ) at 
low x, than that obtained after the standard LO DGLAP evolution. This led to the 
conclusion that the pure LO DGLAP evolution might be not complete at low x. 

On the other hand, the description of polarized data based exclusively on the 
DGLAP evolution has improved after including the DGLAP evolution terms calcu- 
lated at the NLO accuracy []14|,[y|. These predictions described data accurately even 
in the region of low x (x > 10 -4 ). However, the relation between the NLO corrections 
and the ln 2 {l/x) resummation has not been clarified yet. It was expected that at low 
x the NLO corrections contained in part the logarithmic terms but it was not clear to 
what extent those corrections overlapped. 

Here we analyze the interplay between the NLO corrections and the corrections 
obtained from the resummation of double logarithms, ln 2 (l/x). We formulate unified 
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evolution equations for the unintegrated parton distributions, where we include the 
NLO DGLAP terms. Parton distributions are not physical observables and, in general, 
they depend on the factorization scheme applied to the NLO analysis (cf. |5l,|^, |T6|1 ). 
They yield (scheme-independent) physical observables after a convolution with the 
Wilson coefficient functions (cf. Eq. (|5|)). 

In this study we will consider the non-singlet component of the nucleon structure 
function, g± = g\ — g™. This component posseses a property that it does not depend 
on the factorization scheme applied within the family of MS'-like factorization schemes 
TRO. This occurs at the level of parton distributions already. 



We formulate the equations including the logarithmic corrections, ln 2 (l/x), and 
the NLO corrections at the MS factorization scheme. We introduce a unique way 
of performing subtractions in the evolution kernels, so as to avoid double counting 
between the double logarithmic terms and the NLO ones in the overlapping regions 
of the phase-space. Predictions for the non-singlet component of the nucleon spin 
structure function, (x, Q 2 ) and its first moment are then obtained. These results 
are compared with the NLO fit to the experimental data, GRSV'2000, Finally, 
our conclusions are listed. 

2. Unintegrated parton distributions and DGLAP evolution 

Differential equations describing the evolution of standard (integrated) parton dis- 
tributions may be transformed into integral equations, if introducing the unintegrated 
parton distributions. These distributions, fi(x,Q 2 ), are defined as: 

Q 2 dk 2 



A qi (x,Q 2 ) = Aql \x)+J k2 -p- fi(x, k 



where Aqi(x,Q 2 ) denote the integrated parton distributions with Aq^lx) describing 
the contributions coming from the non-perturbative region, Q 2 < k$. The cutoff /cq is 
usually ~ 1 GeV 2 . 

The standard DGLAP equations written for the evolution of the integrated non-singlet 
parton distribution, Aq N $, 

° ;Aq NS (x,Q 2 ) = a s (Q 2 )(AP®Aq NS )(x,Q 2 ), (2) 



dlnQ 2 

with as(Q 2 ) = s ^ - , transform to the following integral equations: 

f NS (x,Q 2 ) = a s (Q 2 )(AP^Aq^ s )(x) + a s (Q 2 ) — (AP ® f NS )(x, k 2 ), (3) 
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if the relation pi) is applied. In particular, we have: — — — AqAx,Q 2 ) = fi(x,Q 2 ) at 

dmQ 2 

the LHS of Eq. (§), and, substituting (g) to the RHS of Eq. @, one obtains the RHS 
of Eqs. (§). 

Coefficients, AP, are the splitting functions calculated for the polarized deep inelas- 
tic scattering. Symbol ® denotes the integral convolution of two functions, (f®g)(x) = 
Jo dyfa dz5(x- yz) f(y) g(z). 

A complete description of the DGLAP evolution of polarized parton distributions 



at the NLO accuracy was performed primarily in one of the MS factorization scheme. 
The splitting functions describing the evolution of the non-singlet components and the 
singlet ones were firstly derived in and they are listed in detail in the Appendix 



of Ref . |n|] . The equations for the DGLAP evolution at the NLO order are structurally 
identical with Eq. (|2|). The only difference is that the non-singlet distribution, Aq^s, 
separates into two parts: the valence one, Aq^ s = Aq — Aq, and the asymmetric one, 
Aq^s — Aq + Aq — (Aq f + Aq'), which have to be evolved with different kernels. Here 
we will evolve the non-singlet isospin-3 parton component, 

Aq NS , 3 = Au + Au - Ad - Ad, (4) 

using the asymmetric MS evolution kernels |1~4| . This component contributes to gi . 



NLO DGLAP evolution of parton distributions depends on the factorization scheme 
applied. In general, scheme dependence dissapears in physical observables, e.g. in 
the structure functions, g\(x,Q 2 ), which are combinations of the parton distributions 
convoluted with the Wilson coefficient functions, AC q ^ g (x): 

9i(x,Q 2 ) = -Y.<{(^ + *s(Q 2 )^C q )®{Aq + Aq))(x,Q 2 ) 

1 q=l 

+ a s (Q 2 )^-(AC g ® Ag)(x, Q 2 )| , (5) 

where 5 is the Dirac 5-function, 5(x — 1). The coefficient functions calculated at the 
MS scheme [|19| are taken from Ref. |14]]- The constant Nf denotes the number of 



active flavours. Here, Nf = 3. 

In case of the non-singlet parton distributions, the dependence on the factorization 
scheme dissapears already at the level of the integrated parton distributions, Aq^ s , 
within the set of the MS'-like schemes [^, [L7|]- The non-singlet component of the 
nucleon structure function is then, obviously, also scheme-independent, and it reads: 

g? s (x,Q 2 ) = \ ({5 + a s {Q 2 )AC q )®Aq N s ;i ){x,Q 2 ) (6) 
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3. Unified evolution equations including ln 2 {l/x) resummation. 



A detailed study of the nucleoli structure function at low x was performed in ||. 
There a complete resummation of the double logarithmic corrections, ln 2 (l/x), was per- 
formed for the unintegrated parton distributions. Those corrections originated from 
summing up the ladder and the non-ladder bremsstrahlung diagrams H,H with the 
requirement of ordering the ratios, k 2 /x n , for exchanged partons. The ratio, k 2 /x n , 
corresponded to the Sudakov variable, (3, the momentum, k n , was the transverse mo- 
mentum, and the fraction, x n , was the longitudinal momentum fraction of the nth 
parton exchanged. The recursive equations for the resummation of the logarithmic 
corrections were then formulated. The integrated parton distributions were obtained 
from the unintegrated ones, using a relation similar to (|l|): 



A qi (x,Q 2 ) = A ft 



(°), 



x 



w 2 dtf 
2 1^ 



fi(x' = x(l + ^),k 2 ), 



(7) 



where the phase-space was extended to W 2 = Q 2 (l/x — 1) corresponding to the total 
energy squared measured in the center-of-mass frame. 

Following 0], we combine the NLO DGLAP evolution and the logarithmic resum- 
mation into a unified evolution equation, including both the DGLAP kernels and the 
ln 2 (l/x) evolution kernels: 

"Q 2 dk 2 



fNs(x,Q ) 



a s (Q 2 )(AP®Aq^ s )(x)+a s (Q 2 ) f 

DGLAP 

- , /Y .. 1 / ' [ Q2/z dk „ 




fNs)(x,k ) 



(8) 



*/*dk? 
(Non — ladder) 

where x' = x(l + ■&). For a detailed form of the kernels see Appendix A. Matrices, 
F§ and Go, represent octet partial waves and colour factors respectively. They are 

described in detail in Appendix A. Symbol 
transform of F 8 /u 2 : 



[z] denotes the inverse Mellin 



F 8 /^ 



5+ioo A,., 

fV"F 8 (o,)/uA 
-ioo 2m 



(9) 
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with the integration contour located to the right of the singularities of the function 
F 8 (u)/u 2 . 

Each of the kernels adds to the homogeneous part of the recursive evolution equa- 
tion. In order to avoid double counting the NLO terms and the ln 2 (l/x) terms in the 
overlapping regions of the phase-space, we propose the following procedure of including 
the NLO terms into the evolution equations (cf. Q). We divide the phase-space of these 
equations into two regions: (i) k 2 < k 2 < Q 2 and (ii) Q 2 < k 2 < Q 2 /z. In the region 
(i) we keep all terms generated by the (non-ladder) double logarithmic corrections and 
add only the regular part of (the LO) and NLO DGLAP terms, AP reg , i.e. this part 
which is not singular at z — > hence does not generate any ln 2 (l/x) contributions. In 
the region (ii) both LO and NLO DGLAP terms do not appear, and we have there 
only contributions from the ladder and the non-ladder double logarithmic terms. This 
procedure is unique, and it uses a proven result of Refs. |12],[R| that the ln 2 {\/x) 



resummation is complete after including the ladder and the non-ladder contributions. 

While integrating the parton distributions over the extended phase-space, k 2 < 
W 2 = Q 2 (l/x — 1) in (0), we may generate singular contributions similar to those 
appearing in the Wilson coefficient, AC q . Therefore, when calculating the non-singlet 
structure function, g± s , we use only the regular part of the Wilson coefficient, 

AC _.i{ a+ ^(!^) + _3_J_ +2+z _g + ^, (1 _ z) }, 10) 

A full form of the Wilson coefficient was taken from Ref. [14||. Symbol () + is defined by 



the following convolution. A function f(z) convoluted with a function (g(z)) + gives: 

1 dz f(z) (g(z)) + = I' dz (f(z) - f(l))g(z) (11) 
o Jo 



4. Results 

We solved numerically the evolution equation (§) for the non-singlet isospin-3 parton 
component, f NS ,3- 

Two different parameterizations of the non-perturbative parton distribution, 

a <1nsA x )' at Qo = 1 GeV2 were used: (0 the GRSV'2000 fit @, and (ii) a simple 
"flat" input 

Apf\x) =Ni{l-x) ni (12) 

with T] Uv = r]d v = 3, r/u = r/s = 7 and r\ g = 5. Normalization constants Ni were 
determined by imposing the LO Bjorken sum rule for Au^ — Ad^\ and requiring 
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that the first moments of all other distributions are the same as those determined from 
the QCD analysis [fZT|. It was checked that parametrization (O) combined with the 



unified equations gave reasonable description of the SMC data on gi S (x, Q 2 ) |§ and 
on gi(x,Q 2 ) [^2|. This fit was also used in [|,|23| to investigate the magnitude of the 



double logarithmic corrections to the spin structure function of proton and to its first 
moment. 

The integrated parton distribution was then obtained by numerical integration of 

fNS,3'- 

Aq NS>3 (x, Q 2 ) = Aq^ 3 (x) + — f NS , 3 (x' = x(l + — ), k 2 ), (13) 



following Eq. (|7|). 



Afterwards, we made a numerical convolution of Aq^s t 3 with the Wilson coefficient 
function, in order to obtain the non-singlet nucleon structure function, gi S , 



g? s (x, Q 2 ) = \({& + a s (Q 2 ) AC g , reg ) <g> Aq NS , 3 ) {x, Q 2 ). (14) 
Figs. H a,b show the results. 



After evolving the GRSV input, the differences between curves resulting from dif- 
ferent approximations (DL+NLO, DL +LO, NLO, LO) were of the order of several 
percent at x > 10~ 3 . The curve representing the unified DL+NLO evolution under- 
estimated the pure NLO results. The same occured also for the evolution of the flat 
input (|i~2|), however, here the difference between the pure LO and other NLO results 
was described by the factor of 2 at x = 10~ 4 . In both cases the differences between the 
DL+NLO, DL+LO and NLO approximations manifested at very low x, x < 10~ 4 . This 
was clearly the region, where ln 2 (l/x) resummation dominated. The results show that 
the shape of the input distribution influenced the evolution of the structure function 
strongly. We checked that the DL+LO and LO results obtained for g^ s are in a good 
agreement with the results of Ref. (see Fig. 5 there). 

Finally, we calculated the first moments of g± s , T^ s , starting from the inputs (i) 
and (ii). We estimated those moments in the following way. We integrated the non- 
singlet structure function, g± s , over x in the region, 10 -4 < x < 1. Since at x < 10~ 4 
the function g± s was not known, we extrapolated g± s into this region with the fit, 
x~ x , and then integrated it over x at < x < 10~ 4 . It is worth mentioning that 
the exponent, A, was found to be 0.33 - 0.44 at Q 2 = 2 - 10 GeV 2 with different 



inputs. That was in agreement with the recent QCD-based estimate from Ref. |24 
A = 0.40 + 0.29. 
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The asymptotic corrections, added to the previous results, changed them at about 
one percent at most. The moments obtained are plotted in Fig. Results obtained 
from the GRSV and the flat inputs, after evolving them with the DL+NLO evolution 
equation, are similar. They follow but underestimate the theoretical predictions for 
T^ s obtained from the QCD analysis at the NNNLO accuracy (Bjorken's sum rule 
||14||). This discrepancy is less then 8 %. In contrast, T^ s estimated from numerical 
integration of the GRSV'2000 fit [|l8j follows Bjorken's sum rule calculated at the 



NLO accuracy, 0(a s ). This clearly suggests that the unified evolution, DL+NLO, goes 
beyond the standard NLO DGLAP analysis, and that the corrections generated by 
double logarithms in the region of low x influence the behaviour of the first moments 
significantly . 

5. Conclusions 

We calculated the non-singlet component of the nucleon structure function, g± s = 
9i ~ 9i-> an d its first moment, using a unified evolution equation. This equation incor- 
porated both the terms describing the NLO DGLAP evolution and the terms which 
contributed to the ln 2 (l/x) resummation. A unique way of including the NLO terms 
into the unified evolution equation was proposed, in order to avoid double counting in 
the overlapping regions of the phase-space. The results obtained from different approx- 
imations were similar, and the minor differences manifested at very low x, x < 10~ 4 , 
where the ln 2 (l/x) resummation dominated. This effect depended on the shape of the 
non-perturbative input distribution, Aq^ s (x). First moments obtained from differ- 
ent non-perturbative inputs evolved with DL+NLO evolution equations were similar. 
These moments followed Bjorken's sum rule calculated at the NNNLO accuracy. The 
discrepancy was at about 8 % at most. This shows that the unified evolution goes be- 
yond the NLO DGLAP analysis. The ln 2 (l/x) corrections influence the Q 2 -dependence 
of the first moments strongly. 



Our results are scheme-independent within the MS'-like family of factorization 
schemes [IT31 (TTI . 
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Appendix A 

Here a brief description of the evolution kernels of Eq. (|8]) is given. DGLAP kernels 
were taken from Ref. fl4| . The full DGLAP kernel AP includes both the LO and NLO 



terms: 

AP = AP LO + a s (Q 2 ) AP nlo . (15) 

r-Q 2 dk 2 

In the homogeneous term, a s (Q 2 ) / -— (AP re9 ® fNs)(x, k 2 ), appearing in (|8|), only 

J k 2 K 

the regular part of the full kernel is included. This is to avoid double counting the 
NLO terms and ln 2 (l/x) terms in the region of the phase-space, k 2 . < k 2 < Q 2 . 

Ladder kernels corresponding to the LO DGLAP kernels taken at z = M generate 
double logarithmic corrections in the region of Q 2 < k 2 < Q 2 / z. 

Non-ladder kernels were obtained in Ref. [|J from the infrared evolution equations 
for the singlet partial waves F , F 8 @,||,[12],[13|]. In H we noticed that extending the 



kernel of the double logarithmic evolution equations from the ladder one, 

a s (Q 2 )AP qq /u, (16) 

to the modified one, 

a s (Q 2 ) (APjco - (F 8 M G )J '(2vrV)) , (17) 

gave a proper anomalous dimension as derived from the infrared evolution equations. 

Matrix Go contained colour factors resulting from attaching the soft gluon to ex- 
ternal legs of the scattering amplitude : 



Gr 







jV 2 -l 
2N 

N 



where iV was the number of colours. 

Further, it was checked that the Born approximation of Fg, 



?^(u)^8n 2 a s (Q 2 )^. (19) 
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gave accurate results for the DL evolution. Martix M 8 was a splitting function matrix 
in colour octet t— channel, 



1_ N F 

2N 2 

N 2N 



(20) 



The inverse Melin transform of F^ orn (uj) then read : 

{z) = 4n 2 a s (Q 2 )M 8 ln 2 (z). 



FBorn 
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(21) 



The evolution equation (§) includes the non-ladder corrections in the Born approx- 
imation (pl|). 
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Figure 1: Non-singlet component of the nucleon structure function, g^ s — g\ — g™, 
plotted as a function of x at the fixed Q 2 = 10 GeV 2 . Thick solid line (DL+NLO) 
shows the results obtained by: (i) numerical solving Eq. (||) for fNS,3 with the input 
given (thick dotted line), (ii) numerical integration of f'NS,3 performed in order to 
obtain Aqws,3, and (iii) numerical convolution of Aq^s,3 with the regular part of the 
Wilson coefficient function (Eq. (|14|)). Thick dashed line (NLO) shows the NLO 
DGLAP evolution of the input. Thin lines correspond to the results: from the ln 2 (l/x) 
resummation including DGLAP terms at the LO accuracy (solid line, DL+LO), and 
from the pure LO DGLAP evolution (dashed line; LO). Upper plot shows the results 
from the GRSV'2000 input ]H| (x > 10 -4 ), lower plot shows the results from the flat 
input (PD (x > 1(T 5 ). 12 



NS 




DL+NLO (GRSV) 
DL+NLO (FLAT) 
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Q 2 [GeV 2 ] 
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Figure 2: First moments of the non-singlet component of the nucleon structure function, 
9i S = 9i - 9i, as a function of Q 2 (2 < Q 2 < 10 GeV 2 ). Thin solid line shows the 
results from the GRSV non-perturbative input, dashed line shows the results from the 
flat input flT2"|). They are compared with the theoretical predictions for the Bjorken 
sum rule (thick solid lines) calculated at the NLO and NNNLO accuracy [ 14 1 , and with 
the first moment estimated from integrating out the NLO GRSV'2000 fit [|T^] (dotted 
line) . 
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